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Abstract. A model for anomalous J/$ suppression in high energy heavy ion collisions is 
presented. As the additional suppression mechanism beyond standard nuclear absorption 
inelastic J /\& scattering with hadronic matter is considered. Hadronic matter is modeled as 
an evolving multi-component gas of point-like non- interacting particles (MCHG). Estimates 
for the sound velocity of the MCHG are given and the equation of state is compared with 
Lattice QCD data in the vicinity of the deconfinement phase transition. The approximate 
cooling pattern caused by longitudinal expansion is presented. It is shown that under these 
conditions the resulting J/^f suppression pattern agrees well with NA38 and NA50 data. 
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1. INTRODUCTION 

The possible existence of a quark-gluon plasma (QGP) phase of highly excited 
hadronic matter is one of the most intriguing questions of high energy physics 
discussed since almost three decades up to now. This novel state of matter has 
been also predicted in lattice QCD calculations (for a review see |l( an d references 
therein) and the critical temperature T c for the ordinary hadronic matter-QGP 
phase transition has been obtained in the range of 150 - 270 MeV depending on 
the number of active quark flavors (this corresponds to the broad range of the 
critical energy densities e c ~ 0.26 — 5.5 GeV/fm 3 ). Since estimates of the NA50 
collaboration for the energy density obtained in the central rapidity region (CRR) 
give the value of 3.5 GeV/fm 3 for the most central Pb-Pb data point, it has been 
argued that the conditions for the existence of the QGP are met in Pb-Pb collisions 
at CERN SPS @- 

The main argument for the QGP creation during Pb-Pb collisions at the CERN 
SPS was the observation of the anomalous suppression of J/\I/ relative yield. J/^J 
suppression as a signal for the QGP formation was originally proposed by Matsui 
and Satz [3j]. The idea is based on the modification of the quark forces in plasma due 
to color screening which would entail a dissociation of c — c bound states, including 
the prominent Jfty at high temperature, observable as a significant depletion of 
the dilepton spectrum around the quarkonium invariant mass. 

The key point of the NA50 Collaboration paper (2j] was the figure, denoted as 



Fig.6 there, where experimental data for Pb-Pb collision values of B ^ aj i^ (the ratio 



<?dy 



of the J/ty to the Drell-Yan production cross-section times the branching ratio of 
the J into a muon pair) were presented together with some conventional predic- 
tions. Here, "conventional" meant that J/\I/ suppression was due to J/^f absorption 
in ordinary hadronic matter. Since all those conventional curves saturated at high 
transverse energy Et, but the experimental data fell from Et — 90 GeV much 
lower and this behavior could be reproduced on the base of J/^> disintegration in 
the QGP 8, this was argued to be circumstantial evidence for QGP formation in 
most central Pb-Pb collisions. In general, the immediate reservation about such 
reasoning was that besides those already known (see e.g. @, @, @, 0, EqI, El] and 



references [12-15] in [2j), other possible models formulated in terms of hadronic 
degrees of freedom could not be ruled out as long as the J/^f absorption rates in 
dense hadronic matter are largely unknown. 



Within a few years the NA50 Collaboration has updated their J/\I/ results [12 
using the nuclear absorption systematics obtained from precise pA results. It ap- 
peared that also these data could be explained with a hadronic comover absorption 
model [13], without the necessity to invoke QGP formation. The question, how- 
ever, arises whether the inputs used in these models are consistent with accessible 
theoretical or experimental information. 

In the present contribution, we present a systematic step towards and general 
description of J/^> absorption in the framework of a statistical analysis which can 
provide a baseline in the search for non-hadronic explanations of anomalous J/\I/ 
suppression. 

The main features of the model introduced in Ref. 15] are 



1. a multi-component non- interacting hadron gas appears in the CRR instead of 
the QGP. All hadrons from the lowest up to Q~ baryon (with their non-zero 
masses) are taken into account as constituents of the matter; 

2. the gas expands longitudinally and transversely; 

3. J/ty suppression is the result of inelastic scattering on constituents of the 
gas and on nucleons of colliding ions. Both "traditional" sources of J/\I/ 
suppression, namely absorption in the nuclear matter and in the hadron gas 
in the CRR, are considered simultaneously. 

The model has the following parameters: the initial time to, the J/\l/-baryon cross 
section cr&, the initial baryon number density n° B , ro in the expression Ra = tqAz 
and the freeze-out temperature T/. 0> . But the last quantity disappears effectively in 
the final estimations because the "natural" freeze-out is enforced by the transverse 
expansion when the rarefaction wave reaches the collision axis. 



2. THE TIMETABLE OF EVENTS IN THE CRR 

For a given A-B collision t = is fixed at the moment of the maximal overlap of the 
nuclei (for more details see e.g. [l3]). As the nuclei pass each other charmonium 
states are produced as the result of gluon fusion. After half of the time the nuclei 



need to cross each other (t ~ 0.5 fm), matter appears in the CRR. It is assumed 
that the matter thermalizes almost immediately and the moment of thermalization, 
to, is estimated to about 1 fm/c [TO, H|. Then the matter begins its expansion 



and cooling and after reaching the freeze-out temperature, T/ ., it ceases as a 
thermodynamical system. The moment when the temperature has decreased to 
T/. . is denoted as tf M ,. Since the matter under consideration is a gas of hadronic 
resonances, no phase transition takes place during cooling. 

For the description of the evolution of the matter, relativistic hydrodynamics 
is employed. The longitudinal component of the solution of the hydrodynamic 
equations (the exact analytic solution for an (1+1) -dimensional case) reads (for 
details see, e.g., 

HQ) ' 



s(t) = , n B {T) = , v g = - (1) 

r t t 

where r = y/t 2 — z 2 is a local proper time, v z is the component of the fluid velocity 
parallel to the collision axis and sq and n° B are the initial densities of the entropy and 
the baryon number, respectively. For n B = and the uniform initial temperature 
distribution with a sharp edge at the border established by the nuclear surfaces, 



the full solution of the (3+l)-dimensional hydrodynamic equations is known [16 
The evolution derived is decomposed into the longitudinal expansion inside a 
slice bordered by the front of the rarefaction wave and the transverse expansion 
which is superimposed outside the wave. Since small but nevertheless non-zero 
baryon number densities are considered here, the above-mentioned description of 
the evolution has to be treated as an assumption in the presented model. The 
rarefaction wave moves radially inward with the sound velocity c s (see Sect. HI). 



3. THE MULTI-COMPONENT HADRON GAS 



For an ideal multicomponent hadron gas in thermal and chemical equilibrium, 
consisting of I species (here, mesons are considered up to K\ and baryons up to 
fi _ ), energy density e, baryon number density n B , strangeness density n s and 
entropy density s are given by (H = c = 1 always) 



1 1 r°° n 2 K- 

Z7T i=i Jo exp^^j^-j+gi 

Zn i=i Jo expj^^- j+gi 
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where = [m 2 + p 2 ) 1 / 2 and m;, Bi, Si, fii, Si and $ are the mass, baryon 
number, strangeness, chemical potential, spin and a statistical factor of species 
i, respectively (an antiparticle is treated as a different species). The chemical 
potential /j, = B^b + Sifis defines that of the overall baryon number, fis and 
that of strangeness, [is- 

To obtain the time dependence of temperature, baryon number and strangeness 
chemical potentials one has to solve the equations (j3J) - (jSJ) numerically with s, tib 
and 77.5 given as time dependent quantities. For s(r) and ub{t) the expressions (CD) 
are taken and rig — since the overall strangeness equals zero during the whole 
evolution, see Sect. [5] 



4. THE SOUND VELOCITY IN THE MCHG 

In the hadron gas the sound velocity squared is given by the standard expression 

c s = ^- (6) 

Since the experimental data for heavy-ion collisions suggest that the baryon number 
density is non-zero in the CRR at AGS and SPS energies 13, 0, Il9], we calculate 



the above derivative for various values of Ub [20j, [21| . 

To estimate initial baryon number density n° B we can use experimental results 
for S-S [13] or Au-Au (l8l . 19| collisions. In the first approximation we can assume 
that the baryon multiplicity per unit rapidity in the CRR is proportional to the 
number of participating nucleons. For S-S collisions we have cLNb /dy = 6 [TtJ and 
64 participating nucleons. For central collisions of lead nuclei we can estimate the 
number of participating nucleons as 2 A = 416, so we have dN B /dy = 39. Having 
taken the initial volume in the CRR equal to ttR\ ■ 1 fm, we arrive at n° B = 0.25 
fm -3 . This is some underestimation because for S-S collisions the beam energy 
was at 200 GeV/nucleon, whereas for Pb-Pb at 158 GeV/nucleon. From the Au- 
Au data extrapolation one can estimate n B = 0.65 fm -3 3]. These values are for 



central collisions. So, we estimate ([6]) for tib = 0.25, 0.65 fm -3 and additionally, to 
investigate the dependence on ns much carefully, for ns = 0.05 fm -3 . The results 
of the numerical evaluation of ([6]) are presented in Fig.[TJ For comparison, we show 
also curves for tib = and for a pure massive pion gas. These curves are taken 
from H. 

The "physical region" we consider between 100 and 200 MeV on the temperature 
axis. This is because the critical temperature for the possible QGP-hadronic matter 
transition is of the order of 200 MeV jlj and the freeze-out temperature should not 



be lower than 100 MeV 18]. At low temperatures we can see a completely different 
behaviour of the cases with ub = and tib 7^ 0. We think that this is caused by 
the fact that for ns 7^ the gas density can not reach zero when T — > 0, whereas 
for ns = it can. For the higher temperatures all curves excluding the pion case 
behave qualitatively in the same way. From T w 70 MeV they decrease to their 
minima (for n B = 0.65 fm" 3 at T = 219.6 MeV, for n B = 0.25 fm~ 3 at T = 202.5 




FIGURE 1. Left panel: Energy density in units of T 4 versus scaled temperature T/T c for 
the present hadron resonance gas model Q compared with the QCD Lattice data [13] and the 
Hagedorn resonance gas with scaled hadron masses to adapt for the case of unphysical quark 
masses in the Lattice simulation [23]. Right panel: Dependence of the sound velocity squared on 
temperature for = 0.65 fm~ 3 (short-dashed), ns = 0.25 fm~ 3 (dashed), = 0.05 fm~ 3 (solid) 
and riB = (long-dashed) . The case of the pure pion gas (long- long-dashed) is also presented. 



MeV, for n B = 0.05 fin" 3 at T = 183.4 MeV and for n B = at T ^ 177.3 MeV) and 
then they increase to merge to each other above T « 250 MeV. 



5. THE COOLING PATTERN VS. SOUND VELOCITY 

In Sect. [3] we have explained how to obtain the time dependence of the temperature 
of the longitudinally expanding hadr on g as. This dependence proved to be very well 



approximated by the expression [20j, [21 



T(r)=T -T- a . (7) 

The above approximation is valid in the temperature range [Ty Q ,T ], where Tf > 
100 MeV, T < T , max and T 0)inax « 230 MeV. We started from T equal to 227.2 
MeV (for n% = 0.65 fm~ 3 ), 229.3 MeV (for n% = 0.25 fm" 3 ) and 229.7 MeV (for 
n° B = 0.05 fm~ 3 ). These values correspond to eo = 5.0 GeV/fm 3 (the initial energy 
density in the CRR has been estimated by NA50 to e = 3.5 GeV/fm 3 ). Then 
we took several decreasing values of T < T 0>max . For every T chosen we repeat the 
procedure of obtaining the approximation ([7]), i.e. the coefficient a of the power 
law. 

We can formulate the following conclusion: in the "physical region" of temper- 
ature and for realistic baryon number densities, the longitudinal expansion given 
by dl|) results in the cooling of the hadron gas described by ([7]) with a = c^(T ), 
namely: 



T(r)=T - M (8) 



where To belongs to the " physical region" . Note that T(r) = T • (^) * is the exact 
expression for a baryonfree gas with the sound velocity constant (for details see 
It should be stressed that c 2 s in (JSJ) is constant depending only on the 
initial temperature To from which the cooling starts. The approximation (jSJ) will 
be used to simplify the evaluation of the J/\& survival factors in the next sections. 

It should be added that for the excluded volume hadron gas model the above- 
mentioned conclusion is no longer valid. We perfprmed the appropriate simulations, 
but approximations of T(r) by a power law turned out to be rather inaccurate and 
yielded a coefficient a which was by a factor two smaller than the squared velocity 
of sound. 



6. J/^> ABSORPTION IN THE EXPANDING MCHG 

As it has been already mentioned in Sect.[2j charmonium states are produced in the 
beginning of the collision, when nuclei overlap. For simplicity, it is assumed that 
the production of cc states takes place at t — 0. To describe J/^ absorption quanti- 
tatively, the idea of Ref. lOj is generalized here to the case of the multi-component 
massive gas. Since in the CRR longitudinal momenta of particles are much lower 
than transverse ones (in the c.m.s. frame of nuclei), the J/^f longitudinal momen- 
tum is put to zero. Additionally, only the plane z = is under consideration. For 
the simplicity of the model, it is assumed that all charmonium states are completely 
formed and can be absorbed by constituents of the surrounding medium from the 
moment of their creation by inelastic scattering through interactions of the type 

cc + h — >D + D + X, (9) 

where h denotes a hadron, D is a charm meson and X stands for a particle which 
is necessary to conserve the charge, baryon number or strangeness. 

Since Pb-Pb collisions are the most relevant case for the problem of QGP search 
(see remarks in Sect. [1]), the further considerations are done for this case. 

It is assumed that the hadron gas, which appears in the space between the nuclei 
after they have crossed each other, also has the shape of the overlap area of the 
colliding nuclei {S e ff) at to m the z = plane. Then, the transverse expansion starts 
as a rarefaction wave moving inward S e ff at to- From the considerations based on 
the relativistic kinetic equation (for details see [El, [lQj]), the survival fraction of 
J/ty in the hadron g clS clS 8b function of the initial energy density eo in the CRR is 
obtained: 

r ( rtfinai Jl^ r cPq Pv<f- ] 

Nh. g .(e ) = J dp T g(p T ,e )-expl- dt 2^J (2^^'^ <Urd,< i£E?J ' ^ 

where the sum is over all species of scatters (hadrons), p v = (E,p T ) and q\ = {E[,q) 
are four momenta of J/^> and hadron species % respectively, cr, stands for the 
absorption cross-section of J/^ — hi scattering, v re i^ is the relative velocity of J/^> 
and hadron h^. M and rrii denote J/^> and hi masses respectively (M = 3097 MeV). 



The function g(pr,^o) is the J/^f initial momentum distribution. It has a gaussian 
form and reflects gluon multiple elastic scattering on nucleons before their fusion 
into a J/^f in the first stage of the collision [25j, [26|, |27|. The upper limit tf ina i of 
the time integration in (fTUi) is the minimal value of (t esc ) and tf, ,. The quantity 
(t esc ) is the average time of the escape of J/^'s from the hadronic medium for 
given values of b and J/^> velocity v = pt/E. Note that the average is taken with 
the weight 

pj/*(r) = T A (r)T B (f-b)/T AB (b) , (11) 

where Tab (6) = J g? 2 s Ta(s)Tb(s — 6), T A (s) = j dzpA(s,z) and p A (s,z) is the nuclear 
matter density distribution taken here as the Woods-Saxon form with parameters 
from [28]. In the integration over hadron momentum in (fTUl) the threshold for the 
reaction (jUJ) is included, i.e. <7j equals zero for {p u + q") 2 < {2m D +m x ) 2 and is 
constant elsewhere (mc is a charm meson mass, mo = 1867 MeV). Also the usual 
Bose-Einstein or Fermi-Dirac distribution for hadron species i is used in (II 01) 



fi{q,t)=f i {q,t) = {2s i + l) exp 



T(t) 



+ 9i 



(12) 



For simplicity, we use (jHJ) as the approximation to T(t) in (|T2l) and and ps{t) 

are solutions of only two equations (j3J) and (jlj) with T given by (jHJ), ns(t) by ([I]) 
and ri5(i) = 0. 

As far as cr, is concerned, there are no data for every particular J — /ij 
scattering. Therefore, we use here universal, energy independent cross sections for 
scattering of charmonia on baryons, cr&, and on mesons, <T m = 2<r b /3, according to 
the quark counting rules. 

As it has been already suggested [7J also J / ^ scattering in nuclear matter should 
be included in any J/^f absorption model. This could be done with the introduction 
of a J/\I/ survival factor in nuclear matter [30] 



[e ) = exp • 



(13) 



where po is the nuclear matter density and L the mean path length of the J/ty 
through the colliding nuclei obtained according to 



PoL{b) 



2T 



AB 



d 2 sT A {s)T B {s-b) 



—T A & + —T B (3-b) 



(14) 



Since the J/ty absorption processes in nuclear matter and in the MCHG are 
separated in time, the J/^f survival factor for a heavy- ion collision with the initial 
energy density eo, could be defined as 



Af(e ) =M n . m {e ) -Nh.g.^o) 



(15) 



Note that since the right hand sides of Eqs. ( fTUl) and (TT3"]) include parts which 
depend on the impact parameter b and the left hand sides are functions of e only, 
the expression converting the first quantity to the second (or reverse) should be 



defined. This is done with the use of the dependence of eo on the transverse energy 
E T extracted from NA50 data (2| (for details see 

To make the model as realistic as possible, one should keep in mind that only 
about 60% of the observed J/\&'s are directly produced during the collision. The 
remainder is the result of \ (~ 30%) and ip' (~ 10%) decays. Therefore the realistic 
J/\I/ survival factor could be expressed as 

7V(eo) = 0.6A/"j/^(e ) +0.3A/" x (e ) + 0.1A/^'(e ) , (16) 



where A/j/^,(e ), J\f x (e ) and A/^'(e ) are given also by Eqs. (TTDT) . (TIB" 
with g(pT,£o), cj/ipN and M changed appropriately (for details see 



and (USD but 
I)- 



7. RESULTS 



Now we can complete the calculations of formula ( JTUl) for values of n° B given in 
Sect. HI The last parameter of the model is the freeze-out temperature. Two values 
T/-. . = 100, 140 MeV are taken here and they agree well with estimates based on 
hadron yields [l8|. The main results are presented in Fig.[2j The original data [2| for 
B ^(J^ b / cr^) b Y b and J/\l/ survival factors given by ( fl~6l) multiplied by B^a^^/ a p ^ Y 
as functions of Et are depicted there (for details see [51]). When comparing the 
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FIGURE 2. J/'f survival factor times B^^a^^/a^y m the longitudinally and transversely 
expanding hadron gas for the Woods-Saxon nuclear matter density distribution and n% = 0.25 
fm~ 3 , 7y o . = 140 MeV, c s = 0.45 and ro = 1.2 fm. The curves correspond to ah = 4 mb (solid) 
and (7b = 5 mb (dashed). The black triangles represent the 1996 NA50 Pb-Pb data, the white 
squares the 1996 analysis with minimum bias, the black points the 1998 analysis with minimum 
bias 0, and the white stars to the 2000 NA50 Pb-Pb reanalyzed data [re- 
present model with the data from the NA50 collaboration in Fig. [21 we want to focus 
on the region below Et ~ 100 GeV, which corresponds to initial energy densities 
below e ~ 3.0 GeV/fm 3 (T < 210 MeV), where the MCHG gives an acceptable 
description of the equation of state for hot and dense matter, in agreement with 
lattice QCD data, see the left panel of Fig. [TJ In our phenomenological analysis 
shown in Fig. [2J, we assume universal cross sections for mesons and baryons, with 



the appropriate thresholds for their dissociation reactions but energy-independent, 
0& = 4 mb (solid line) and 5 mb (dashed line). Note that the solid line perfectly 
describes the data between the onset of the deviation from the nuclear absorption 
baseline at Et ~ 40 GeV and the limit of applicability of the MCHG picture at 
E T rsj 100 GeV. Does this result disprove the claim of the NA50 collaboration that 
the onset of anomalous suppression at Et ~ 40 GeV is evidence for the creation 
of a new form of matter, made of deconfmed quarks and gluons (QGP) in this 
experiment? 

Such a claim would be premature. The weak point is the question about the 
charmonium dissociation cross sections. Although the absorption cross section 
for J/ip in cold nuclear matter has been measured in pA collisions by the NA50 



experiment to be ctj/^jv = 4.2 ±0.5 mb [3JJ, this may be a result of the downfeeding 



from Xc and ip' with larger cross sections, so that the "true" absorption cross 
section for the ground state component is between 2 and 3 mb, see also (32|. No 
experimental information for the charmonia absorption cross sections on higher 
baryonic resonances or on mesons is available yet. 

However, there is progress in theoretical approaches. Based on a diagrammatic 
approach to quark exchange processes in hadron-hadron scattering [33| . the disso- 
ciation cross section for charmonia on mesons have been calculated [34J to have a 
sharp rise at threshold to maxima between 1 ad 10 mb, depending on the channel, 
followed by a fast decrease in energy due to vanishing overlap integrals between 
asymptotic mesonic states. A recent calculation within a fully relativistic approach 
has confirmed the result for J/ip dissociation by pion impact [35]. Similar calcula- 
tions for baryon impact [36] show that both assumptions made in this work could 
be disproved: the universality and energy independence. 



8. CONCLUSIONS 

Properties of the hadron resonance gas model under conditions reached in the 
ultrarelativistic heavy ion experiments are very close to the phase transition region 
estimated from the lattice simulation data, see Fig.[IJ We have shown that with J/ty 
absorption cross sections on hadrons of 4 mb, an overall satisfactory description of 
NA38/NA50 data on J/ty suppression could be given here. Therefore, one could 
argue that the NA50 Pb-Pb data do not provide evidence for the production of 
deconfined matter in the central rapidity region of a Pb-Pb collision. 

It is an interesting task for forthcoming work to use those energy dependent, 
non-universal cross sections of the quark exchange model in the calculation of 
charmonium dissociation in the MCHG model. We expect that the outcome of 
such a calculation will leave room for the discussion of in-medium modification 
(increase) of dissociation rates by changes in the spectrum (broadening) of the 
final state open charm hadrons which determine (effectively lower) the reaction 
thresholds. Examples have been given for the ir — p gas [37| and for nuclear matter 



381 ] . A heuristic extension of the MCHG to include spectral broadening due to 



the Mott dissociation of hadrons (Mott-Hagedorn resonance gas) has been given 



m 



391 ] . but deserves a microscopic foundation. Therefore, we conclude that the 



hadronic absorption cross sections of the J/^f need to be determined to a higher 
accuracy before the anomalous J/^f suppression could be interpreted as a good 
signal for QGP formation in central heavy-ion collisions. 
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